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We analyze the possibility of observing new physics effects in the b — » d penguin 

■ 

■ amplitudes. For this purpose, we consider the decay mode — ► K°K°, which has 
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only b — > tZ penguin contributions. Using the QCD factorization approach, we find 
very tiny CP violating effects in the standard model for this process. Furthermore, we 
show that the minimal super symmetric standard model with LR mass insertion and 

Qh! R-parity violating super symmetric model can provide substantial CP violation effects. 

Oh. Observation of sizable CP violation in this mode would be a clear signal of new physics. 
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1 Introduction 



Recently there have been a lot of interests to look for new physics effect beyond the standard 
model (SM). The recent measurement of the indirect CP violating parameter S^x s in the 
decay mode B® — > (pKs, which is a pure b — > sss penguin induced process, may provide the 
first indication of new physics beyond the SM [1]. Within the SM, the mixing induced CP 
asymmetry in the B® —> 4>K S mode is expected to be equal to that of B® — > tpK s [2] within a 
correction of 0(\ 2 ). The most recent updated data on S^ko by BABAR [3] agrees within one 
standard deviation with the value of S^k s whereas, the Belle data [4] has about 2a deviation. 
Therefore, the presence of new physics (NP) in this mode has not yet been ruled out from 
the available data. In principle, the new physics can affect either the B% — -B° mixing or the 
decay amplitudes. Since the new physics effect in the mixing can affect equally to both the 
cases the above deviation may be attributed to the decay amplitude of B® — > (f)K s or more 
generally to the b — > s penguin amplitudes. The next obvious question is: Do the b — > d 
penguin amplitudes also have significant new physics contribution. The present data does 
not provide any conclusive answer to it. The obvious example is the B® — > nn processes, 
which receive contribution from b u tree and from b — > d penguin diagrams. The charge 
averaged branching ratios of all the three processes 5° — > 7r + 7r~, B% — > 7r°7r° and B ± — > 7r ± 7r° 
[5] and the CP violating parameters C n7T and S nn in B% — > it + it~ process [6] have already 
been measured. The present situation is: the measured branching ratio for the color allowed 
process B® — > 7r + 7r~ is about two times smaller than the QCD factorization calculation and 
the measured B® — > 7r°7r° color suppressed branching ratio is about six times larger than 
the corresponding QCD factorization calculations [7]. Thus the discrepancy between the 
theoretical and the measured quantities imply the following two consequences. 

• The QCD factorization may not be a very successful theory for the charmless B decays. 

• There may also be significant new physics effect in the b — > d penguins as speculated 
in b — > s penguins. 

Recently Buras et al. [8] have shown that the observed B — > nn data can be explained 
in the standard model if one includes the large nonfactorizable contributions. In this paper 
we would like to look into the second possibility i.e., the existence of any new physics in 
b — > d penguin amplitudes and indeed if it does, could it be detectable. For this purpose, 
we consider the decay mode B% — > K°K° which has only b — > d penguin contribution. The 
significance of this decay mode is that it originates from b — > dss penguins with dominant 
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contributions coming from the QCD penguins. If one assumes that the penguin topology is 
dominated by internal top quark, the CP asymmetry parameters would vanish in the SM. 
However, as pointed out by Fleischer [9], contribution from penguins with internal up and 
charm quark exchanges are expected to yield nonzero CP asymmetry in — > K°K° mode. 
Thus, the study of CP asymmetries in this mode may provide an interesting testing ground 
to explore new physics effects. The CP averaged branching ratio has recently been measured 
by the BABAR collaboration [10] 



which agrees with the SM predictions [7] . Although, the measured branching ratio does not 
provide any indication for a possible new physics effect, the measurements of CP violation 
parameters in near future will certainly establish/rule out the possibility of new physics in 
the b — > d penguin amplitudes. This decay mode has recently been analyzed by Fleischer and 
Recksiegel [11] in the SM. They have shown that this channel may be characterized through 
a surface in the observable space from which one can extract the relevant information. In this 
paper we consider the impact of new physics effect on the CP violation parameters. We show 
that the minimal supersymmetric model with LR mass insertion and the supersymmetric 
model with R-parity violation can provide significant CP violation effect, the observation of 
which would be a clear signal of new physics. 

The paper is organized as follows. In section 2, we discuss the basic formalism of CP vio- 
lation. Using QCD factorization approach, we estimate the CP averaged branching ratio and 
CP violating parameters in the SM. The effects of new physics on the CP violating parameters 
are discussed in section 3. The contributions arising from minimal supersymmetric standard 
model with mass insertion approximation and from R-parity violating supersymmetric model 
are discussed in subsections 3.1 and 3.2 respectively. Our conclusion is presented in section 
4. 

2 CP violating parameters in the Standard Model 

We first present a very general treatment of the CP violating parameters. The time depen- 
dent CP asymmetry for B® — > K°K° can be described by 




(1) 



a KK (t) = C KK cos AM Bd t + S KK sin AM Bd t , 



(2) 
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where 

1-|A| 2 2Im(A) 

Ckk - YTW' KK ~ ~T+W' { ) 

In the above expression A corresponds to 

g A(gj -> K°K°) 

p A(B d — > K°K°) { ) 

where, g and p are the mixing parameters defined as 



Q 



M* — -T* 

M 12 2 1 12 (5) 



p \| m 12 - |r 12 ' 



The off diagonal element of the mass matrix is given by the matrix element of the AB = 2 
transition as 

(B° d \H AB=2 \B° d ) = M 12 - l -Y l2 . (6) 

In the standard model, the box diagrams are dominated by the H^-boson and top quark in 
the loop, as a result of which, one obtains (ignoring terms of 0(Ti 2 /Mi 2 )) 

q _ VgVtd _ 2i/3 (7) 

p v tb v; d ~ e ■ [n 

The amplitude for the decay mode B° d — > K°K°, which receives dominant contribution in 
the SM from QCD penguins with top quark in the loop can be written as 

A(B° d - K°K°) = V tb V t * d P t , (8) 

where are the CKM matrix elements which provide the weak phase information and P t is 
the penguin amplitude arising from the matrix elements of the four quark operators of the 
effective Hamiltonian. The amplitude for the corresponding CP conjugate process is given 
as 

A(B° d - K°K°) = V t * b V td P t . (9) 

Thus one gets 
and hence 

Ckk = S kk = 0. (11) 

Thus if the measured CP violating asymmetries in B° — > K°K° deviates significantly from 
zero it would be a clear signal of new physics. However, the decay amplitude also receives 
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some contribution from the internal up and charm quarks in the loop. Therefore, the CP 
violating parameters may not be zero identically. Now including the effects of u, c, t quarks 
in the loop and using CKM unitarity one can write the decay amplitude as 



MB, 



K°K°) 



^uPut + \Pct — ^uPut 



(12) 



where \ q = V qb V* d , P qt = P q — Pt are the penguin amplitudes, 5 = 5 ct — S ut is the relative 
strong phase between them and 7 is the weak phase. The parameter r is defined as 

Pa 



1 

Rb 



ut 



with R b 









H; 


j 





pi _|_ ^2 



(13) 



Thus one obtains the CP asymmetries as 

sin 2a + 2r cos 5 sm(2(3 + 7) — r 2 sin 2/3 



Srk — 
Crk = 



I _|_ r 2 _ 2 r C os 5 cos 7 
— 2r sin 5 sin 7 



(14) 



1 + r 2 — 2r cos 5 cos 7 

where a, /3 and 7 are the three angles of the unitarity triangle. Thus, to know the precise 
value of the CP violating asymmetries we must know the values of r and 5. The CP averaged 
branching ratio for the process is given as 



K U K U ) + Br(5^ -> K°K ' 



where the individual branching ratios are given as 

TB°\Pc.m\ 



87rmR 



A(B° d -> K°i?°) 



(15) 



(16) 



We now use the QCD factorization approach to calculate the branching ratio and CP asym- 
metry parameters. The effective Hamiltonian describing the process under consideration 
is 



n 



eff 



10 



(17) 



£ cp 3 + c g o g 

-i=3 

where q — u, c. O3, • • • ,Oq and O7, • • • , Oiq are the standard model QCD and electroweak 
penguin operators respectively, and O g is the gluonic magnetic penguin operator. The values 
of the Wilson coefficients at the scale ji m m b in the NDR scheme are given in Ref. [12] as 

C 3 = 0.014 , C 4 = -0.035 , C 5 = 0.009 , C 6 = -0.041 , C 7 = -0.002a 



C s = 0.054a , C 9 = -1.292a , C w = 0.263a , C g 



-0.143 . 
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In the QCD factorization approach [7], the hadronic matrix elements can be represented in 
the form 



0r>0| 



l + 5> n c£ + C(AqcdM) 



(19) 



where (K K |Oj|-B^)f act denotes the naive factorization result and Aqcd ~ 225 MeV is the 
strong interaction scale. The second and third terms in the square bracket represent higher 
order a s and Aqco/^b corrections to the hadronic matrix elements. 

In the heavy quark limit the decay amplitude for the B d — > K°K° process, arising from 
the penguin diagrams is given as 



G f 



V £ q=u,c 



q Q 10 , q , ( q a l . q 

a 4 ^~ T a 10a + r x I °6 ^" + a 8a 



X , (20) 



where X is the factorized matrix element. Using the form factors and decay constants defined 
as [13] 



(K (p K )\sYb\B° d (PB)} 



(K°(q)\dYl 5 s\0) 



(Pb+PkT ~ 



q 



(21) 



we obtain 



X = (K (^)|s7 M (l-75)6| J B0(p B ))(K (g)|rf 7 At (l-75) S |0) 
= -ifi<Fo(m 2 K ) (m| - m 2 K ) . 



(22) 



The coefficients af s which contain next to leading order (NLO) and hard scattering correc- 
tions are given as [14, 15] 



C 4 + -^ + 



C 3 a s C F 



N 4vr N 



C 3 [F K + G K (s d )+G K (s b )] 



+ C 1 G K (s q ) + (C A + C 6 )'£G K (s f )+C g G Kt! 

f=u 

C5 a s C F 
N 4tt N 



a\ = C 6 + — + — 



C 3 [G' K {s d ) + G' K {s b )} + C l G' K {s q ) 



+ (C 4 + C 6 ) E G' K (s f ) + C g G' Kfj 

f=u 
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at = C 8 + 



9l 
3 ' 



a 8a = 



a s C F 
4tt N 



(Cs + C w )\ £ e f G' K (s f ) + «| [e^( Srf ) + e b G' K {s b )\ 

f=u 



l \0a 



Ct s C F 

4vlv 



3 b 3 
(C 8 + Cio)- E e / G ^( s /) + ^2 [ e ^( s ^) + e b^(s fe )] 



(23) 



where q takes the values u and c, N — 3, is the number of colors, CV = (iV 2 — 1)/2N. The 
internal quark mass in the penguin diagrams enters as s/ = to 2 /to|. The other parameters 
in (23) are given as 

Fx = _i 2 ] I1 A_i 8 + /£ + /£ > 

/■i 1 — 2a; 

f K = / dx g(x)<f) K (x) , g( x ) = 3- lnx-3i7r, 

jo 1 — x 

f K = —T^K77w^J -M*)J -M*)J -My), 



N F B ^ K (0)m 2 B 



G 



K,g 



X 



>K 



Ga-(s) = ^-in-^ + 4 ( l dx<p K {x) 
6 6 rrih Jo 



3 3 nib 

x / c^u w(l — u) In [s — u(l — u)(l — x) — ie] , 
Jo 



rl 3 3 

G^ >9 = -y dx-(fP K {x) = -- , 



- ln-^- + 3 / f/.c 
to 6 



X 



/ w(l — u) In [s — m(1 — u)(1 — x) — ie] 
Jo 



(24) 



The light cone distribution amplitudes (LCDA's) at twist two order are given as 



i B (x) = N B x 2 (l - x) 2 exp - 



,2 „2 



2ul 



(25) 



0^(x) = 6x(l - x) , ^(x) = 1 > 

where iVs is the normalization factor satisfying Jq 1 dx<j)B(x) = 1 and uj b = 0.4 GeV. The 
quark masses appear in G(s) are pole masses and we have used the following values (in 
GeV) in our analysis 



TO,, 



TOrf = to s = 0, m c = 1.4 TO5 = 4.8. 



r x = 2m 2 K /{m b — m s ){m s + ma) denotes the chiral enhancement factor. It should be noted 
that the quark masses in the chiral enhancement factor are running quark masses and we 
have used their values at the b quark mass scale as m b (m b )=AA GeV, m s (m b )=90 MeV and 
m d (m b ) = 6.4 MeV. 

For numerical evaluation we have used the following input parameters. The value of 
the form factor at zero recoil is taken as F (0) = 0.38, and its value at q 2 = m 2 K can be 
obtained using simple pole dominance ansatz [13] as F (m 2 K ) = 0.383. The values of the 
decay constants are as fx = 0.16 GeV and fs = 0.19 GeV, the particle masses and the 
lifetime of B® meson t b o = 1.536 ps are taken from [5]. Thus we obtain the amplitude (in 
units of 10~ 2 ) 



G F 



A(B°^K°K°) = 



'71 



(26) 



V ub V: d (13.56 + % 4.59) + V cb V c * d (14.98 + i 2.06) 

V ub V: d (14.32 e m - 7 °) + V cb V: d (l5.12 e^ 83 °)" 

We use the values of the CKM matrix elements at the la CL in the Wolfenstein parameter- 
ization from Ref. [16] as 

A = 0.226518:8^, A — 0.80118:820 > P = 0.18918Z, ^ = 0.35818:82, (27) 

which correspond to the angles of the CKM unitarity triangle 

sin 2a = -0.1418:2 , P = (23.8^)° , 7 = (62lg)° ■ (28) 

With these input parameters we obtain the CP averaged branching ratio as 

B(B° d -> K°K°) = (9.15 ± 0.30) x 10~ 7 , (29) 

which is slightly below the central experimental value (1). Since, we have used the LCDA's 
at twist two order our predicted result is slightly lower than that of Ref. [7] where they have 
included the twist three power corrections in the distribution amplitudes. From Eqs. (12) 
and (26) one can obtain 

r w 2.6 and S w 11° . (30) 
With these values we get the CP asymmetry parameters in the SM as 

S KK = 0.061 C KK = -0.163 . (31) 
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By allowing the CKM matrix elements to vary within their la range as given in Eqs. (27) 
and (28), we obtain the correlation between Skk and Ckk in the SM as shown in Figure-1, 
which gives the constraints 

0.02 <S KK < 0.13, - 0.17 < C KK < -0.15 . (32) 

Thus, if the measured values of CP asymmetry parameters will be outside the above ranges, 
would be a clear sign of new physics. 

-0.15 
-0.152 
-0.154 
-0.156 
-0.158 

o - 016 

-0.162 
-0.164 
-0.166 
-0.168 
-0.17 

0.02 0.04 0.06 0.08 0.1 0.12 0.14 

Figure 1: The correlation plot between Skk and Ckk for the B Q d — > K°K° process in the 
SM, where we have used r = 2.6, 5 = 11° and the CKM parameters are varied within the 
range as given in Eqs. (27) and (28). 

3 New Physics effects on the CP violating parameters 

Here, we consider the effect of new physics on the CP violating parameters. Because of 
the new physics contributions, the CP asymmetry parameters (14) become modified. In 
principle, the new physics can affect either the B® — B d mixing or the decay amplitudes. Let 
us first investigate its effect in the mixing phenomena. In the presence of new physics, there 
are additional contributions to the mixing parameters arising from the new box diagrams. 
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These contributions to the AB = 2 transitions are often parametrized as [17] 



1/12 =r m e*», (33) 



V Mf 2 M 

where M\i = M-j| M + is the off diagonal element of the mass matrix, contains contri- 

bution both from the SM and from new physics. Hence, the ratio of the mixing parameters 
q/p as given in Eq.(5) becomes 

V 

Thus, in the presence of new physics, the mixing induced CP asymmetry in B® — > ipK s can 
be given as 

S^ Ks = sin(2/3 + 20 m ) (35) 

However, since the present world average on the measurement of S^k s = 0.726 ± 0.037 [18] 
agrees quite well with the SM prediction S$k s = 0.715^o;q45 [19], we do not consider the 
effect of NP in mixing in our analysis. 

Now we consider the effect of new physics on the CP violating parameters arising from 
the new contribution to the standard model decay amplitude (12). In the presence of new 
physics the amplitude can be written as 

A(B° d - K°K°) = A SM + A NP = \ u P ut [l - re^ + r NP e w »e^} , (36) 

where tnp = \A^p/ X u P ut \, On and 5n are the relative weak and strong phase between them. 
The amplitude for the corresponding CP conjugate process can be obtained by changing the 
sign of the weak phases. 

Thus the CP asymmetry parameters (3) become 



where 



Sr K = ~y, and C££ = -f , (37) 



X = sin 2a + 2r cos<5 sin(2/5 + 7) — r 2 sin 2/3 + 2tnp cos 5n sin(2a; + On) 
—2rr NP cos(<5 — 5n) sin(0 N — (2/3 + 7)) + r NP sin(2o; + 20 n) , 

Y = 1 + r 2 + r NP — 2r cos 5 cos 7 + 2r NP cos 5 N cos N 
—2rr NP cos(5 — 5n) cos(7 — On) , 



Z = 2 



r sin S sin 7 — tnp sin Sn sin On + ttnp sin(<5 — 5n) sin(7 — On) 



(38) 
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The branching ratio in the presence of new physics is given as 



Br(B° d -> K°K°) = Br SM ( 1 + 



A 



NP 



i-SM 



+ 2r 



NP 



A 



NP 



i5M 



COS(/> 



N 



(39) 



where Br SM denotes the SM branching ratio and 0jv is the relative phase between the new 
physics and standard model amplitudes. 

Now, we consider two beyond the standard model scenarios: the minimal supersymmetric 
standard model with mass insertion approximation and R-parity violating supersymmetric 
model and study their effects on CP violation parameters in the following subsections. 



3.1 Contribution from minimal supersymmetric standard model 
with mass insertion approximation 

Here, we analyze the decay process B® — > K°K°, in the minimal supersymmetric standard 
model (MSSM) with mass insertion approximation. This decay mode receives supersymmet- 
ric (SUSY) contributions mainly from penguin and box diagrams containing gluino-squark, 
chargino-squark and charged Higgs-top loops. Here, we consider only the gluino contribu- 
tions, because the chargino and charged Higgs loops are expected to be suppressed by the 
small electroweak gauge couplings. However, the gluino mediated FCNC contributions are of 
the order of strong interaction strength, which may exceed the existing limits. Therefore, it 
is customary to rotate the effects, so that the FCNC effects occur in the squark propagators 
rather than in couplings and to parameterize them in terms of dimensionless parameters. 
Here we work in the usual mass insertion approximation [20, 21] where the flavor mixing 
% — > j in the down-type squarks associated with g# and (Ja are parametrized by (5^ B )ij, with 
A, B = L, R and i,j as the generation indices. More explicitly (Sf L )ij = (V^ M^V^)ij/rri~, 
where Mj is the squared down squark mass matrix and trig is the average squark mass. 
is the matrix which diagonalizes the down-type quark mass matrix. 

The new effective AB = 1 Hamiltonian relevant for the B — > K°K° process arising from 
new penguin/box diagrams with gluino-squark in the loops is given as 



n eff = ~^/2 U 



(40) 



E {cro, + era) + c? p o a + c? p d s 

.i=3 

where (O g ) are the QCD (magnetic) penguin operators and Cf p , C^ p are the new Wilson 
coefficients. The operators Oi are obtained from Oj by exchanging L <-> R. 
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To evaluate the amplitude in the MSSM, we have to first determine the Wilson coefficients 
at the b quark mass scale. At the leading order in mass insertion approximation the new 
Wilson coefficients corresponding to each of the operator at the scale /i ~ m ~ M w are 
given as [21, 22] 



cr 



AG F V tb V* d m\ ^ LL 



C 



NP 
4 



AG F V tb V* d mj 



riNP 

°6 



AG F V tb V* d m\ 



(«0 

(«0 



C, 



NP 



AG F V tb V* d m\ 
2 v / 2vra s 



13 



13 



13 



-^(x) - b -B 2 {x) - lp x (x) - \p 2 {x) 
-\b x (x) + \b 2 {x) + ±P x (x) + ^P 2 (x) 
^B 1 (x) + ±B 2 (x) - -Lp^x) - l -P 2 {x) 
+ ^P 2 (x) + ip^x) + ^P 2 (x) 



2G F V tb V td m\ 



t d LL ) 13 {^M,{x) - l -M 4 {x) 



X) 



(41) 



where x = m|/m~. The loop functions appear in these expressions can be found in Ref. [21]. 
The corresponding Ci are obtained from Cf p by interchanging L <-> R. It should be noted 
that the (5^) 13 contribution is enhanced by (mg/m b ) compared to the contribution from 
the SM and the LL insertion due to the chirality flip from the internal gluino propagator in 
the loop. 

The Wilson coefficients at low energy C^ p (fi ~ m b ), can be obtained from Cf p {Mw) 
by using the Renormalization Group (RG) equation as discussed in Ref. [12], as 



C( f i) = U 5 ( l x,M w )C(M w ) , 



(42) 



where C is the 6x1 column vector of the Wilson coefficients and U 5 (/z, M w ) is the five-flavor 
6x6 evolution matrix. In the next-to- leading order (NLO), U 5 (/i, M w ) is given by 



u 5( „, m w) = (1 + iMj\ u«v m w) { 1 - *mi 



An 



An 



(43) 



where \jf\n, M w ) is the leading order (LO) evolution matrix and J denotes the NLO 
corrections to the evolution. The explicit forms of U 5 (/i, M w ) and J are given in Ref. [12]. 

Since the O g contribution to the matrix element is a s order suppressed, we consider only 
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leading order RG effects for the coefficient Cg P ', which is given as [23] 

r *P(m ) - ( \ 2/21 ( as ( mt A 2/23 (A A) 

For the numerical analysis, we fix the SUSY parameter as rriq = m g = 500 GeV, a s (Myy) = 
0.119, a s {m b = AA GeV)=0.221, a s {m t = 175 GeV)=0.107. Now substituting the values of 
the RG evoluted Wilson coefficients C/^^m^'s in Eq. (23) we obtain the corresponding a^s 
and hence with Eq. (20) the amplitude. Assuming that all the mass insertion parameters 
(Sa B ) 13 have a common weak phase, we obtain the fraction of new physics amplitude as 
defined in Eq. (36) 

r NP ~ 0.33 (|(^ L ) 13 | - \{5 d RR )i3\) +465.86 (|fe) 13 | - |(4l)i 3 |) ■ (45) 

It should be noted that because of the opposite chiral structure of the currents O; and O h the 
LL and RR and also the LR and RL contributions occur with opposite sign. As seen from 
Eq. (45), the LR(RL) insertions have dominant effect because of the m g /mb enhancement. 
We use the limits on the (5f L )i3 and (<^r)i3 mixing parameters from [24] for x — 1 as 

|(#l)is|<0.2 K^DisI < 0.01 (46) 

and assume that only one of these gives a dominant SUSY contribution. This gives (tnp) ll < 
6.6 x 10~ 2 and (tnp)lr < 4.66. Since the new physics effect due to LL insertion is almost 
negligible it will not provide any significant effect on the CP violating observables. The 
correlation plot between and for LL insertion is shown in Figure-2, where we use 
r = 2.6, 5 = 11° as obtained from QCD factorization analysis, the central values of the 
CKM weak phases from (28), the relative weak phase 9 N = n and vary the relative strong 
phase 5^ between and 2n. In this case, because of the negligible new physics contribution 
one gets only tiny CP violating effects. In Figure-3, we present the correlation plot for LR 
insertion, where we have used (tnp)lr = 4.66, < 5n < 2tt and a representative set of weak 
phases 9 N = ir, tt/2, 7r/3, 7r/4. For r and 5, we have used the values as obtained from QCD 
factorization (30). As expected, in this case large CP violation can be generated. 

The branching ratio (39) versus 4>n is plotted in Figure-4 for \A np /A S m\ = 2.19. One 
can see from Figure-4 that the observed data can be easily accommodated in minimal su- 
persymmetric standard model with LR mass insertion. 

Thus in future, if sizable CP violation effects will be observed in — > K°K° mode, the 
minimal supersymmetric standard model with LR mass insertion may be a strong contender 
of new physics to explain the data. 
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Figure 2: The correlation plot between Sg% and Cg% for the B° d -> K°K° process in the 
MSSM with only LL insertion, where we have used (r NP ) LL = 6.6 x 1CT 2 , the weak phase 
9 N = 7r, r = 2.6, 5 = 11° and varied the strong phase S N between and 2n. 

3.2 R-parity violating supersymmetric contribution 

We now analyze the decay mode in the minimal supersymmetric model with R-parity vi- 
olation (RPV). In the supersymmetric models there may be interactions which violate the 
baryon number B and the lepton number L generically. The simultaneous presence of both 
L and B number violating operators induce rapid proton decay, which may contradict strict 
experimental bound. In order to keep the proton lifetime within the experimental limit, one 
needs to impose additional symmetry beyond the SM gauge symmetry to force the unwanted 
baryon and lepton number violating interactions to vanish. In most cases this has been done 
by imposing an ad hoc symmetry called R-parity defined as, R = (— i^( 3B+L+2S \ where S 
is the intrinsic spin of the particles. Thus R-parity can be used to distinguish the particles 
(R = +1) from their superpartners (R— —1). The conservation of R-parity implies that the 
supersymmetric particles must be produced in pairs and the lightest supersymmetric particle 
(LSP) must be stable. However, there is no compelling reason to require the conservation of 
R-parity. Less restrictive symmetries- conservation of baryon/lepton number alone can be 
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Figure 3: The correlation plot between Sg% and Cg% for the B° d -> K°K° process in the 
MSSM with LR mass insertion, where we have used r = 2.6, 5 = 11°, (r NP ) LR = 4.66, a set 
of weak phases 9n = n, ir/2, n/3, 7r/4, and varied the strong phase 5^ between and 2n. 



imposed to prohibit the unwanted proton decay. Extensive studies has been done to look for 
the direct as well as indirect evidence of R-parity violation from different processes and to 
put constraints on various R-parity violating couplings [25]. 

Here, we consider only the lepton number violating effects. The most general -R-parity 
and lepton number violating super-potential is given as 



Wy - -XijkLiLjEl + \' ijk LiQjDl , 



(47) 



where, i, j, k are generation indices, Li and Qj are SU (2) doublet lepton and quark superfields 
and E%, D k are lepton and down type quark singlet superfields. 

Thus the relevant four fermion interaction induced by the R-parity and lepton number 
violating model is 

1 



2m? V 



-8/A) 



+ K32Ki2(dal»Lbp)(s p ^Rs a ) + \* 3 \' i21 (d a ^Rbp)(sp^Ls a ) 



(48) 
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Figure 4: The branching ratio of B® — ► K°K° process (in units of 10~ 6 ) versus the phase 
0iv (in degree). The horizontal solid line represents the experimental central value and the 
dotted lines represent the la range. 

where 77 = " S ^~P an( j fl Q — \\_l n ^ The QCD correction factor r]~ 8 /P° arises due to running 
from the sfermion mass scale mj, (100 GeV assumed) down to the m b scale. 

The amplitude for B° d — > K°K° process in the RPV model is given as 



A n (B° - K°K°) = W 8//3 °^ 

ny ' 8ml 



^ (y y* — y y* \ 

\ A i22 A H3 A i31 A 122 J 

2 m K ( \i \l* _ w \l* \ 

™ V™ 1 ™ \ l A i21 A i23 A i32 A H2j 

[m b - m s ){m d + m s ) v 7 



(49) 



where we have kept only the leading order factorization contributions. We use the parame- 
terization \' i2 2 A 'ii3 = - A i3i A '*22 = ke ' e and A- 32 A'*i2 = -A- 2 iA'* 23 = he 10 , assuming the same 
weak phase for all the RPV couplings. The limits on the couplings |A- 32 A'*i2l = |A- 21 A'*23l 
are obtained from B® — > 07r decay in Ref. [26] 

h = \Ks2 A '*n2\ = \ A 'i2i A m\ < 4.0 x 10- 4 . (50) 
In our analysis we use k — k± < 4.0 x 10~ 4 and obtain the new physics parameter 

r NP < 3.92 . (51) 
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The correlation plot between C^k an d $kk f° r the above value of r NP is shown in Figure-5, 
for some representative values of the weak phase and < 5 NP < 2ir. The values of r and 
5 are used as derived from QCD factorization. In this case also one can get observable CP 
violation effects. Plotting the branching ratio (39) vs. 4>n for \A np /A S m\ = 1-84 we can see 
from Figure-4 that the observed branching ratio can be easily accommodated in the RPV 
model. 




4 Conclusion 



The recent measurement of the mixing induced CP asymmetry in — > <\>Ks which has 
significant deviation from sin(2/3)^K- g may provide the first indication of new physics effects 
present in the b — > s penguin amplitudes. In this paper, we have investigated the possibility 
of observing new physics effects in the b — > d penguin amplitudes. We have considered 
the decay mode B% — > K°K° which proceeds through the quark level FCNC transition 
b — > dss, receiving contributions only from one-loop b — >• d penguin diagrams. If one would 
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assume only the top quark exchange in the penguin loop as usually done, the CP asymmetry 
parameters would vanish in the SM. However, contributions from penguins with internal 
up and charm quark exchanges are expected to yield small non-vanishing CP asymmetries. 
Thus, if significant CP asymmetries will be found in this channel then it would be a clear 
indication of new physics effects in b — > d penguin amplitudes. However, as discussed 
in Ref. [27], the nonfactorizable long-distance charm penguins may also give significant 
contributions which in turn yield sizable CP asymmetries. In that case it is practically 
impossible to disentangle the new physics effects from the nonfactorizable charm and GIM 
penguins without any additional assumptions. However, very recently, it has been pointed 
out by Beneke et al [28] that the nonfactorizable charm penguin contributions are of higher 
order in l/m b expansion. Thus the observation of sizable CP asymmetry in this mode may 
be considered as the signal of new physics. 

Using QCD factorization approach, we found the CP averaged branching ratios in the 
SM for B® — > K°K° process as ~ 0.9 • 10~ 6 , which is slightly below the present experimental 
value. The CP asymmetry parameters are found to be Skk — 0.06 and Ckk — —0.16. 
Allowing the CKM parameters to vary within their la limits, we obtained the allowed ranges 
as 0.02 < Skk < 0.13 and —0.17 < Ckk < —0.15. If the observed values would deviate 
significantly from the above ranges would be a clear signal of new physics. We next analyzed 
the decay mode in the MSSM with mass insertion approximation and found that the LR 
insertion has significant effects than the LL or RR insertions. In this case one can have 
significant CP violating asymmetries. Considering the R-parity violating supersymmetric 
model we found that one can also obtain significant CP violation with the present available 
RPV couplings. Therefore, the future experimental data on B^ — > K°K° CP violating 
parameters will serve as a very good hunting ground for the existence of new physics beyond 
the SM and also support/rule out some of the existing new physics models. 
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